Josephson tunnelling spectroscopy of negative U centers 
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We consider a superconductor-insulator-superconductor (SIS) junction in which the tunnelling 
through the insulating barrier is dominated by a localized "negative U" center. We show that 
the LcR product of the junction depends sensitively on the spectrum of impurity states, and in 
near resonant condiditions exhibits an anomalously large IcR product which can exceed the famous 
Ambegoakar-Baratoff limit by an arbitrarily large factor. The analysis is extended to problems 
in which there is an array of negative U centers in the junction. We also discuss general reasons 
to expect significant violations of the optical conductivity sum rule in most SIS junctions and of 
the Ambegoakar-Baratoff result when the superconductors emerge from a non-Fermi liquid normal 
state. 



It has long been known that there exist localized defect 
states in solids that can be characterized as "negative U 
centers." What this means is that the effective interac- 
tions are such that, 



U = E{2) + Ei0) - 2Eil) < 



(1) 



where E{n) is the energy of the state when occupied by 
n electrons, i.e. one can view these impurities as "pair 
binding" regions. The existence of such centers has been 
found to explain many of the anomalous properties of var- 
ious materials in which an impurity or defect with an odd 
number of electrons produces a diamagnetic (rather than 
the naively expected paramagnetic) respoBse. In particu- 
lar, this model was discussed by Andersonl in connection 
with dangling bonds in chalcogonide glasses. 

An obvious question arises whether the effective at- 
traction between electrons implied by Eq. (1) can be, in 
some way, harnessed to provide a mechanism for pairing 
in a superconductor. The trouble, of course, is that the 
negative U centers with which we are familiar are highly 
localized, and the materials in which they occur are good 
insulators rather than supcrconducto*K. However, in a 
multicomponent system, it is possiblea for electron itin- 
eracy to derive from one component which is proximity 
coupled to a set of localized negative U centers, from 
which the pairing arises. Indepjd, superconductivity with 
Tc — 1.5K has been observecB in semiconducting PbTe 
doped with around 1% Tl, a welLknowno negative U cen- 
ter. Along these lines, two of ua3 have recently proposed 
that negative U centers between the copper-oxide planes 
may play a role in enhancing Tc in certain cuprate high 
temperature superconductors, most notably in the Hg 
based materials. 

In the present paper, we address the less ambitious, 
but related issues: What is the effect of weak coupling 
between a conventional superconductor and a negative 
U center? Can Josephson tunnelling be used as a spec- 
troscopic probe to detect the presence of negative U cen- 
ters? Manifestly, in the weak coupling limit, the presence 
of a negative U center can have little or no effect on the 



strength of the superconducting state, itself. However, as 
we will show below, the incipient pairing on the negative 
U center leads to an anomalous enhancement of the IcR 
prodjHct; it can exceed the famous Ambegoakar-Baratoff 
limitO by an arbitrarily large factor. By sweeping the 
chemical potential on the impurity site (perhaps by ap- 
plying a suitable gate voltage) the impurity level can be 
moved closer and further from resonance, which leads 
to predictable spectrocopic variations of the critical cur- 
rent. We also discuss the generalization of these results 
and their possible pertinence to recent experiments in 
layered superconductors. 

The Case of a single negative U center: To be- 
gin, we consider a tunnel junction between two pieces of 
bulk superconductor through a barrier containing a sin- 
gle negative U center. Specifically, we consider the model 
Hamiltonian 



H — Hl + Hji + Hj, 



(2) 



where Hl and Hn are the Hamiltonians of the left 
and right superconductors, which we take to be iden- 
tical superconductors and well |-|approximated by the 
Bougoliubov-deGennes cquationsQ, i.e. 

Hl = - ^^)LlL,, + HL^Llk-a + h.c). (3) 



and the negative U center is described by 

Hu = -|i7|cjc|cjcx + {e- y)^cUa- (4) 

a 

Here, c\ creates an electron on the impurity site with 
spin polarization a and and L\ create, respectively, 
an electron at the edge of the junction in the right and 
left superconductors. The tunneling Hamiltonian (which 
we treat as a small perturbation) is 



Htun = - 'Yy^LL^Ca + IrR^C^ -f h.C.]. 



(5) 
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(Of course, other interactions may be important under 
some circumstances, including direct tunnelling through 
the junction (not involving the impurity) and more com- 
plicated pair-tunelling and other interaction terms; under 
appropriate circumstances, the present model includes 
the most important interactions, and it is sufficient for 
our present purposes.) 

In the absence of the tunnelling term the eigenstates of 
the impurity are singly occupied, unoccupied or doubly 
occupied (in a spin singlet). The ground state is always 
one of the latter two, so it is more convenient to express 
the results in terms of the single and two-particle exci- 
tation energies ei and 62, both positive, rather than to 
work directly with U and s. For example, for e — > 
and \U\ < 2(e — ^), the groundstate has no particles in it 
and £1 = e — /i and £2 = 2{e — fx) — \U\. Although all of 
our results can be derived for the entire range of param- 
eters we are specifically interested in the case where ei is 
large. In the single impurity case we will use ei/A ^ 1 
to simplify the expressions. 

It is now a straightforward exercise to compute the zero 
temperature Josephson coupling to lowest (4i/i) order in 
powers of Hum'- 



J = A 



TrNptLtB. 



Alog'(2ei/A) 



£2 



(6) 



where Np is the normaL state density of states at the 
Fermi energy. The signu of the J is such as to favor 
alignment of the phases {0l and Or) across the junction; 
I.e. the Josephson energy is Ej — — Jcos(0l — Or)- As 
usual, the critical current is simply /j, — 2eJ/ch. 

The Josephson coupling is the sum of two terms: The 
first term is the contribution to the pair tunnelling across 
the junction from processes in which first one electron 
and then the other passes from one side of the junction 
to the other. The second term comes from processes in 
which the second electron enters the junction before the 
first has left it. In this sense, one can think of the first 
term as being the tcxtbookEl result with an effective tun- 
nelling interaction Tlr = titfl/ei. However, the sec- 
ond term, which involves correlated pair tunnelling, is a 
qualitatively new effect. Being proportional to A^ it di- 
verges at the point at which the unoccupied and doubly 
occupied impurity states are degenerate (the impurity 
becomes partially occupied), i.e. a resonance condition 
for pairhopping. Note that, because the junction con- 
ductance is proportional to the square of the amplitude 
for single-particle transmission across the barrier, there is 
no analagous correlated pair-tunnelling term which con- 
tributes to 1/R. 

By itself, this expression is not terribly illuminat- 
ing. However, as usual, we can express J in units of 
the normal-state junction conductance, 1/R, so that the 
highly detail dependent factors of the tunnelling matrix 
elements cancel, and we are left with a quantity, with 
units of energy, which in some sense measures how big 
the Josephson coupling is in natural units. R can be com- 



puted using a straightforward generalization of standard 
perturbative methods. It is independent of how far the 
impurity is from resonance (£2). In the limit of small ^, 
the result, the principal result of this paper, is 



elcR 



ttA 
~2~ 



1 



A 



log2(2ei/A) 



(7) 



Note that far from resonance, when the second term is 
absent, this result reproduces the standard Ambegoakar- 
Baratoflo relation. However, near the resonance, where 
the unoccupied and doubly occupied states of the impu- 
rity are nearly degenerate with each other, the negative 
U center gives rise to a dramatically enhanced IcR prod- 
uct. When the resonance condition is too nearly satified, 
i.e. e2/A < Npt^/ci, the perturbative expression is no 
longer useful. 

Luckily, this single impurity problem can be solved for 
arbitrary £2 (but still perturbatively in tunnelling): in- 
tegrating (projecting) out the singly occupied impurity 
states and the two superconductors reduces the problem 
to a rather simple effective Hamiltonian 



(8) 



for the unoccupied and doubly occupied states in the ef- 
fective proximity field 

^=M£M^(tie-^-|-4e-«) (9) 

and e'^P is the contribution from the first term in Eq. ^ 
to the Josephson energy. The ground state energy 



Ej{0l~Or) 



£2 



(10) 



can now be used not only to obtain the altogether larger, 
0{t'^), contribution right at the resonance (£2 = 0) 



- - ^^-l°g(2£r/A) /,4 + ,4^ + 2(t,t,)^cos(^.-^^«), 



E:r = 



ei 



(11) 



but also to smoothly connect with Eq. ^. 

In principle, when this kind of an effectively single im- 
purity junction is realized, the characteristic energies, £2 
and £1, can be varied by applying a voltage in the junc- 
tion, so that the resonant condition can be tuned. In 
many circumstances one is faced with a collection of such 
impurities and so it is worth generalizing the derivation 
for an extended system. 

Tunnelling through a correlated region: The 
general structure of the perturbative calculation is un- 
changed. For the single impurity the result is determined 
by purely local correlations (both of the superconductors 
and, trivially, the impurity) and therefore a single param- 
eter, — , controls the physics. For an extended barrier, a 
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far more diverse roster of possibilities exists due to the 
interplay of disorder and coherence scales of the super- 
conductors and in the barrier. 

In the tunnelling Hamiltonian we now allow for a non- 
uniform (though still short-ranged) hopping: 

Htun = f dr[tL{7^Li{r)c,{r) + tR{r)Ri{r)cA,^ + h.c.]. 

(12) 

The superconductors are still described via Bogoliubov- 
deGennes Hamiltonian or, more generally, by their 
anomalous (single particle) green functions !FL/ji{r,T)u. 



The impurity region is described entirely by its 
(time-ordered) two particle correlation function 
(r^cx(l)c|(2)c|(3)c|(4)), where 1...4 is a short-hand 
for space and time coordinates, (fi,ri). 

We now make a crucial simplifying assumption in de- 
scribing the nonsuperconducting layer. We will take its 
single electron green function to be short ranged in space 
and time. This is certainly well justified when there is 
a gap or, more generally, whenever the relevant electron 
energy scales are large compared to A. 

More specifically, we can always define the pair (P2) 
and partidle-particle propagators (Pi and Pq) by decom- 
posing the two particle correlation function as 



(r,q(l)cT(2)c (3)c (4)) = G(l, 2)G(3, 4)^2(1, 2; 3, 4) + G(l, 4)G(2, 3)Pi(l, 4; 2, 3) + G(l, 3)G(2, 4)Po(l, 3; 2, 4) (13) 



where G(ri,r2,Ti — T2) = {TrC^-^{rT)c]{f ',r')) is the 
imaginary time-ordered Green function, and where Pi — *■ 
0, P3 0, and Pq ^ 1 at large distances. The pertru- 
bative expression for the Josephson coupling involves in- 
tegrals over this quanatity. We will simplify this expres- 
sion by making two physically motivated assumptions. 
Firstly, we will assume that the integral of G and any 
other less strongly peaked function / can be simplified 
as follows 



dT' / dr'G{r,f',T')f{f',T') 



1 1- 
-f{r,-, 
e e 



(14) 



where e is a characteristic single particle excitation en- 
ergy (analogous to ei above). If /(r, 0) is finite, we 
can let /(r, 1/e) /(r, 0), but if / is weakly singu- 
lar as r 0, 1/e plays the role of a cutoff - see be- 
low. Among other things, this permits us to substitute 
P„(l,2;3,4) ^ P„(l,l;4,4) = P„(l,4) in Eq. |l|. Sec- 
ondly, since we are interested in a barrier region which 
supports significant pairing fluctuations, we will assume 
that all the interesting correlation effects are reflected in 
the behavior of P2, an will consequently make the sim- 
plifying approximations Pq = 1 (no significant charge 
density wave fluctuations) and Pi = (no significant 
magnetic fluctuations). 

With these simplifying approximations, we can bring 
the analog of Eq. |^, J — Ji + J2, into a somewhat more 
manageable form with 



Ji = 2 j drdr'TLR{7^TLR{r') j\T\T{f-f',T)\^ 



2|^(o,i; 



(15) 



drdf'tl{r)tl{r') / drPi (f, f', r), 




where TLR{r) = tL{r)tR{f)/e. Notice that the first term, 
as before, only depends on the parameters of the tunnel- 
ing region through the typical single particle energy, e, 
and, so can be viewed as a direct tunnelling term but with 
an effective hopping matrix, Ti^R(f). More importantly, 
the hopping elements, tL{r) and tR{r)^ enter Ji and J2 
very differently, and so disorder will generally have a very 
different effect on the two contributions. 

To make our analysis concrete, we adopt a simple 
model in which there are random variations of (r) and 
tj^(r), but no correlations between and tR: 



tL(r)ti(r') = ti[aL + (1 - az.) exp(-|f - r'|/a)] (17) 



(and similarly for Ir), but tL{f)tR{r') = 0, where <L(r) = 
tL\/oiL signifies the configuration average. For simplicity, 
we ignore all other sources of disorder. 

In terms of this model, it is possible to obtain expres- 
sions for the various contributions to the Josephson cou- 
pling in terms of various averages over the various prop- 
agators. These expressions are simple in terms of the 
implicitly defined intrinsic coherence lengths 



drd''"VPi(f,0,T) - ^2 



d-l 



£2 

drd'^-V|.F(f,0,T)|2 = ^!l^[AiV;.]2, 



(18) 
(19) 



(16) 



where, in addition 

dTT{Q,T) = and ^'^(O, ^) = ATV^ log[e/A]. 

(20) 

It is interesting to remark that the various expressions in- 
volving !F, which we have evaluated here in the context 



3 



of BCS mean-field theory, depend on the behavior of the 
superconducting leads at distances short compared to ^o- 
Here, as we discuss below, the correlation functions inti- 
mately reflect the fact that the normal state of the leads 
is a Fermi liquid. 

We do not display the complete result, here, since it 
is long. In the limit of weak disorder, a — *■ 1, the con- 
figuration averaged Josephson energy per unit area (area 
= L''^^) is readily evaluated 



J 



2^(^NFtLtRf 



Co 



d-1 



log2(6/A) . (21) 



Here, the ratio of the normal contribution and correlated 
pair tunnelling contributionsv to the Josephson coupling 
is proportional to (e2/A)(^o/C2)'^^^; unless the pair cor- 
relations are at quite low energy and extend over quite 
long distances, the anomalous ( J2) pair-tunnelling term 
will be insignificant. 

Now we consider a large disorder limit, a — + 0, and for 
simplicity assume as well that £_l, <C ^2, ^o- In this 
case. 



J 



2cLCRA{TrNFtLtR) 



Ti 



(22) 



+T2^V^log2(6/A) 



where Tj are geometric factors of order 1. Here, so long 
as ^2 ^ £,L, £,R, there is a large enhancement of the corre- 
lated pair tunnelling term. The physics of this enhance- 
ment is very simple - if ^2 is large, a pair can tunnel into 
the barrier where tR is large, then propagate along the 
barrier, and finally tunnel out where is large. This 
spatial structure amplifies the sort of resonance effects 
we found in the single impurity case. 

Discussion: With regards to application of these 
ideas to layered superconductors, our conclusions are in- 
triguing, but not unambiguous. Large IcR products have 
recently been measurpi along the c-axis in various lay- 
ered superconductorscj, as have dramatic shifts of op- 
tical spectral weight over large ranges of energies. The 
two phenomena are thought by many to be related. We 
find thatp-all simple models of SIS junctions (and more 
generallyO) exhibit large spectral weight shifts leading 
to apparent sum rule violations as the system enters the 
superconducting state. 

To see this, consider the case of a simple junction, i.e. 
Htun — — '^„[tL'lR<j + h.c.]. Rather than computing the 
actual optical conductivity we use the single-band ver- 
sion of the optical sum rule (which relates the integrated 
spectral weight in a given band to the expectation value 
of the kinetic energy) to compute the integrated oscillator 
strength for the normal junction and the superconduct- 
ing junction. The result is that 



dLj[aN{Lj) - as{uj)] = 4|i|2iV|.Alog 



W 
A"' 



where and as a-re, respectively, the frequency de- 
pendent junction conductance with superconducting and 
non-superconducting leads. Since the f-sum rule must 
ultimately be satisfied, presumably the correct interpre- 
tation of this result is that spectral weight is transferred 
over large energy ranges, beyond those described by our 
model, i.e. on energy scales on the order of the band- 
width, which is obviously much larger than A. Since this 
occurs even in the absence of an IcR anomaly, we are 
at present unable to make a clear statement concerning 
the expected spectral consequences of such an anomaly. 
Thus, the observed apparent sum-rule violations cannot 
be definitively linked with the presence or absence of in- 
teresting correlations in the barrier region. And, unfor- 
tunately, the experimental data on large IcR products in 
the high temperature superconductors, while very strik- 
ing indeed, may potentially be hard to interpret as well, 
as it is intertwined with pseudo-gap phenomena which 
onset well above the superconducting transition temper- 
ature. 

Perhaps yet a potentially more vexing concern has to 
do with the non-universality of the IcR product even 
when tunnelling through an uncorrelated insulator. The 
problem is that the result is strongly influenced by the 
short distance/time properties of the superconductor. In 
weak coupling these are inherited from the "normal" 
state. We have already seen one manifestation of this: 
the logarithm in the anomalous term is inherited from the 
underlying Fermi Liquid (this is the same log responsi- 
ble for the superconducting instability). More generally, 
even the simple proportionality, J- ^ ANp, very familiar 
from the BCS theory (and in part responsible for the Am- 
begoakar/Baratoff result), between the order pacameter 
(JT) and the mean field (A) need not be generial3 either 
deep enough in the superconducting phase or when the 
normal phase is not Fermi Liquid like. 

In summary, we have considered a problem of tun- 
nelling through a single negative U impurity and an ex- 
tended region with incipient pairing. We find, under 
favorable circumstances, anomalous enhancement of the 
Josephson coupling due to correlated pair tunnelling pro- 
cesses. 
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